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Abstract. The scattering of an electromagnetic plane wave by a spherical 
particle was solved analytically by Gustav Mie in 1908. The Mie solution 
is expressed as a series with very many terms thus obscuring the physical 
interpretations of the results. The purpose of the paper is to try to illustrate 
this phenomenon within the Ramsauer framework used in atomic and nuclear 
physics. We show that although the approximations are numerous, the Ramsauer 
analytical formulae describe fairly well the differential and the total cross sections. 
This allows us to propose an explanation for the origin of the different structures 
in the total cross section. 
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1. Introduction 



The subject of light scattering by small particles is present in several scientific areas 
such as astronomy, meteorology or biology [1]. The first model by Lord Rayleigh in 
1871 dealt with light scattering by particles whose dimensions are small compared 
to the wavelength. A significant improvement came with the Mie solution [2] 
which describes light scattering by spherical particles of any size. This extension 
is important in astronomy and meteorology where light can go through aerosols which 
are particles dispersed in the atmosphere. In Section [21 we study the Mie series for 
the light scattering on dielectric spheres in air. In an apparently different domain, 
the scattering of a low energy electron by atoms studied by Ramsauer showed 
surprising structures, and it took several years before the solution was imagined by 
Bohr: describe the electron as a plane wave [ll[5l[6]! The Ramsauer effect is certainly 
the first phenomenon showing the wave properties of matter. This framework has since 
then been used to describe very different types of collisions such as atom-atom [7] or 
even neutron- nucleus [HI UHl IHl El] • Our idea is that since light behaves like a 
wave, it should be also possible to apply Ramsauer's ideas to the light scattering by 
dielectric droplets. The Ramsauer effect is described in Section[3l Finally, in SectionlH 
we compare the predictions of Mie and Ramsauer for the total scattering cross section 
of light over a rain drop. 



2. The Mie predictions for light over a sphere of non absorbing dielectric 

The Mie solution is detailed in the [Appendix A| Let R be the radius of the sphere, 
let n be the index of refraction of the dielectric, and let crtot be the total cross section. 
In the context of light scattering, the extinction efficiency factor Qo = crtot /{t^R^) is 
often used. 

Figure [1] (a) shows, for three values of the index of refraction, Qc as a function 
of the size parameter x = 2'kR/X = /cout-R: where fcout is the wave number of light 
outside the sphere. Each curve is characterized by a succession of maxima and of 
minima with superimposed ripples. The amplitude of the large oscillations and of the 
ripples grows with n. In the next sections we will show that all these features can be 
understood with the Ramsauer approach. 



3. Description of the Ramsauer effect 



The Ramsauer effect was discovered in 1921 [5] while studying electron scattering over 
Argon atoms. The total cross section versus the electron energy showed a surprising 
dip around 1 eV. In Figure [2] we see recent measurements of electron over Krypton 
and neutron over Lead nucleus. 

The idea to model the phenomenon is to consider the incident particle as a plane 
wave with one part going through the target and another which is not (Figure [3]). The 
two parts recombine behind the target and then interfere with each other, producing 
the oscillating behaviour. Depending upon the impact parameter b, light rays going 
through the drop accumulate a phase shift (Figure[4|). The calculation of Qo.r is given 
in the [Appendix B[ 
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Figure 1. The extinction efficiency factor versus (a) x and (b) Ar = 2x (n — 1), 
for non absorbing spherical particles with relative refractive indices n = 1.05, n = 
1.33, and n = 2.00. x is given by the relation x = 2itR/X = koutR- Ramsauer 
solution for n = 1.05 is also given in low panel. In (a) and (b), the vertical scale 
applies only to the lowest curve, the others being successively shifted upward by 
2. 



where Ar = 2R{kin — fcout) = 2R{n — 1) fcout = (71 — 1). Note that the first 
multiphcative factor of Equation [T] explains the fact that the amphtude of the large 
oscillations grows with n. The extension efficiency factors, when plotted against Ar, 
show a universal shape (see Figure [TJ(b)). Our model and the Mie prediction are in 
good agreement for the three refractive indices used in Figure [T] 

A small fraction of light, internally reflected (IR) twice, will also contribute to the 
forward flux. The amplitude of the internal reflection coefficient r = (n — l)/(n + 1) 
is small (r = 1/7 for visible light upon water drop in air). The phase shift will be 
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Figure 2. (a) Electron-Krypton normalized total cross section versus energy [6]. 
(b) Normalized Neutron- Lead total scattering cross section versus energy |12j . 



roughly = 2 x 2 R x fc;,! — ARnkout — ^nx. Thus the main Ramsauer extinction 
efficiency factor Qo.r is modulated by the internal reflection factor 

^IR = 1 + cos Air = 1 + cos(4 n x). (2) 

Finally, the global extinction efficiency factor is 

Qc = Qc,R X i^iR. (3) 
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Figure 3. Qualitative picture of the Ramsauer plienomenon. The wavelength 
of the light is supposed to be reduced in the dielectric. This picture shows the 
case where the contraction of the wavelength between (a) outside and (b) inside 
the medium (dashed line) is such that they come out in phase. Thus the sphere 
becomes invisible resulting in an almost zero cross section. 
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Figure 4. Definition of the variables used in the text to calculate the extinction 
efficiency factor of light over a spherical particle of radius R. 

4. Application of the Ramsauer approach to the scattering of Ught over a 
drop of water and comparison with the Mie solution 

If we consider visible light ray (A = 0.6 /.jm) on a drop of water (n = 1.33) in air, 
the Mie and the Ramsauer predictions are compared in Figure [5] The Mie solution 
is a series with more than 600 terms that has to be computed for each value of 
the abscissa. On the contrary, the Ramsauer approach is a purely analytical function 
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(Equations ([T|). ^) with clear physical concepts. The Ramsauer model reproduces 
quite well both the amplitude and the peak positions of the Mie prediction, except 
for the main peak position which is lower by almost 10%. At small parameter 
values, the two curves differ. The analytical formula (U), obtained under the light 
ray approximation, is not expected to be a good description of the reality when the 
wavelength of the incident light is very much larger than the size of the droplet. This 
is the Rayleigh regime where the cross section behaves as l/A**. In terms of the 
variable x, this implies, as x 1, a. x'^ behaviour whereas our formula approaches a 
parabola. Even though the small ripples are also present, they look more attenuated 
in the Ramsauer curve. Another particularity of our model is that it justifies the ratio 
between the Ramsauer-pseudo period and internal reflection-pseudo period. Since one 
has already determined the phase difference for each case, it is straightforward to 
derive their ratio using the equations that require constructive interferences for the 
fcout i?-axis, 

^ = (4) 
Air 2n ^ ' 

In the case of a raindrop with n = 4/3, there is thus a factor 8 between the two 
pseudo periods. Therefore, we have a simple physical explanation for the origin of the 
oscillations at two frequencies, and we have derived the ratio of their periods. 




AF, = 2x(n-1) 



Figure 5. Comparison of the extinction efficiencies, the Mie prediction (thin 
line) (calculated according [f 31 ). and the prediction of the Ramsauer model and 
multiple internal reflections (thick lino), x is given by the relation x = 2ttR/X = 

^out 

The Ramsauer framework allows also, through the Huygens-Fresncl principle, to 
calculate the unpolarized differential cross section. This is developped in [Appendix C| 
where Ramsauer and Mie are compared on a particular example. 

5. Conclusion 

We have shown that within the experimental errors the Ramsauer effect predicts an 
extinction efficiency factor comparable to the one given by the Mie solution. The large 
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oscillatory behaviour is understood as the consequence of the interference between the 
fraction of light going through and the fraction avoiding the drop. The origin of the 
small ripples can be traced back to the internally reflected light interfering with the 
other two components. 
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Appendix A. The Mie solution 

The theory describing the scattering of an electromagnetic plane wave by a 
homogeneous sphere was originally presented by Gustav Mie |2] . Particles with a size 
comparable to the wavelength of visible light are relatively common in nature. The 
laws describing the total scattered intensity as a function of the incident wavelength 
and the characteristics of the particle are much more complex than the one for the 
Rayleigh scattering. The 1/A^ dependence of the total scattered intensity in the 
Rayleigh case is not true anymore in the general case for the Mie solution. Thus the 
extinction cfhciency factor Qc = crtot/ {ttR'^) depends on the radius and the relative 
refractive index of the particle. 

The Mie theory uses Maxwell's equations to obtain a wave propagation equation 
for the electromagnetic radiation in a three dimensional space, with appropriate 
boundary conditions at the surface of the sphere. The extinction efhciency factor 
obtained is 

2 °° 

Qe = — ^(2^+l)Rc(af + 6,). (A.l) 

The Mie scattering coefficients ag and hg are functions of the size parameter 
X = 2'kR/X = koutR and of the relative index of refraction n, 

^ xMx)i''giy) - yi^'i{x)My) 2) 
<eixKi{y) - yCiix)Ceiy) 

^ yMx)i^'g{y) - x'ip'g{x)My) 3^ 
yQ{x)C'i{y) - xCg{x)Ce{y) 

where y = nx. ipi{z), Qi^) are the Riccati-Besscl functions (the prime denotes 
differentiation with respect to the argument) related to the spherical Bessel functions 
ji{z) and yg{z) through the equations 

Mz) = zjg{z), (A.4) 

O(^) = z.ie{z) - ^zyg{z). (A. 5) 

Numerically, an infinite sum cannot be computed. Thus it is necessary to truncate 
the series and keep enough terms to obtain a sufficiently accurate approximation. The 
criterion developed by Bohrcn CF in [15j was obtained by extensive computations. 
The number of required terms N has to be at least the closest integer to x + 'ix'^^^ + 2. 
For instance, for a raindrop of 50 fim radius and a visible wavelength of 0.6 /xm, the 
number of required terms is = 558. Nowadays, the computers have reached a point 
where the computing time is no longer a problem, nevertheless it is interesting to 
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compare to an approximate closed form method in order to obtain some insights of 
the different physical processes. The plots of the Mie solutions are obtained using 
the Fast Mie Algorithm of Pawel Gliwa [T3] for total cross sections and the MiePlot 
program of Philip Laven [Tl] for differential cross sections. 



Appendix B. Determination of the extinction efficiency factor: the optical 
theorem 



In the scattering theory, the wave function far away from the scattering region must 
have the form 



ikr 



f{0)- 



The optical theorem reads as 
47r 

fTtot = — Im/(0 = O), 



(B.l) 



(B.2) 



where f{6 = 0) is the forward scattering amplitude. Under the approximation for 
the scattering of a scalar (spinless) wave on a spherical and symmetric potential, the 
scattering amplitude at a given polar angle 9 can be written as a sum over partial 
waves amplitudes, each of different angular momentum £ as follow 



oo 



(B.3) 



1=0 



where rji is the inelasticity factor (77^ = 1 in our case of a non absorbing sphere), 5^ is 
the phase shift {5i is real for a pure elastic scattering), and Pi{cos6) is the Legendre 
polynomial. 

From Figure[4l 5i = (fcjn — kont)Rcos^. Under the approximation of the forward 
scattering, P^(cos6') 1. Let us introduce the impact parameter h = Rsmip, such 
that £ = fofcout (from the Bohr momentum quantization : bp = £h and p = hk). By 
substituting the discrete sum over ^ by a continuous sum over £ or over the impact 
parameters 



E 



d£ 



db, 



the forward scattering amplitude can be written in term of the impact parameter 



/(0 = 0) = — 
i 



^i2{kin — kaut) R cos 4' 



- 1 bdb, 



or in term of the angle tp 



f{d = 0) = 



2i 



_ ei2(fc.„-feo„t)flcos^j j^,Qg2 ^_ 



(B.4) 



(B.5) 



Then the expression for the total cross section is obtained from Equation (|B.4p 
and setting w = cos tp, 



Utot = Im 
= Im 



i47ri?2 



w 



1 - cKfc.n-fco„t)2flu,\ 



i 2ttR^ - i AnR^ 



,i(fc„ 



)2Rw 



dw 



(B.6) 
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Then, from integration by parts, we obtain 



ftot 



Im 



Im 



i27ri?2 - i47ri?2 



i 27ri?^ - i inR 



gi(A;i„-fcout)2-R-iii 

W- 

i(^in ^ fcout) 2i? 

gi(fci„-fco„t)2_R 







gi(fci„ — /Cout)2ii-u) 

i(fcin - kout) 2-R 
fco„t)2i?. _ 2 



i(/Sin - fcout) 2i? 

sin(2i?(fcin - fcout)) 1 



(2i?)2 (fci„ - A^out)^ 



(B.7) 



2_R (fcin fcout) 



Sin(i? (fcin - fcout)) 



If we introduce the parameter Ar 
the cross section becomes 



2 

2i?(fci 



(Jtot = 27ri?^ 



1 



sin Ar / sin Ar/2 



At 



Ar/2 



i? ( fc'in fcout ) 

~ fcout), the analytic expression for 

(B.8) 



In this expression, the total cross section approaches 2ttR^ at high energies. It 
is the so called "black disk" approximation. Nevertheless, from a purely geometrical 
viewpoint, a collision occurs if the distance between the two sphere centers is less 
than ri + r2, where ri and r2 are the radii of the two spheres. So the geometrical 



cross section is equal to 7r(ri 



r2 



In our case, the photon may be considered as 



a particle with a diameter equal to its reduced wavelength A = A/27r. Consequently, 
Equation (jB.Sp becomes 



(Ttot = 27r i? 



A 

47r 



1 



sin Ar / sin Ar/2 
Ar/2 



(B.9) 



Finally, an expression for Ar as a function of the index of refraction n is needed. 
According to Maxwell approach, it is fc'in = n fcout- With this relation and the fact 
that A = 27r/fcout 



(Ttot = 27ri?^ 1 



n- 1 



1 - 2 



sinAR / sin Ar /2 
Ar/2 



or 



(B.IO) 



Qe,R = 21 



Ar 



1 - 2 



sIuAr /sin Ar/2 



Ap 



Ar/2 



This result is derived under the approximation of a scalar light whereas the Mie 
solution has to do with the full vectorial light. Therefore our formula applies only to 
unpolarized light scattering. 



Appendix C. Ramsauer solution for the differential cross section 

Let us first have a reminder about the Huygens-Fresnel principle. The propagation 
of light can be described with the help of virtual secondary sources. Every point of 
a chosen wavefront becomes such a secondary source and the light amplitude at any 
chosen location is the integral of all these sources emitting spherical waves towards 
that observation point. The secondary sources are driven by the incident light. Their 
amplitude is 

1 1 + cos^ 

2 ^incidcntQ'-' 
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where A is the wavelength and ^incident is the amphtudc of the incoming hght. 9 is 
the angle between the direction of the incident light and the direction from the virtual 
source to the observation point. At last the elementary area where the secondary 
source stands is dS. 

We decide to use the plane (P), shown in Figurc|4l as the location of our secondary 
sources. The shadow of the water sphere on the plane (P) is the disk CD. If the plane 
(Q) is chosen as the origin of the phases, the incident light is a plane wave which, 
on the plane (P), has the value e'^'^^. The amplitude distribution of the secondary 
sources can be split into two parts: an undisturbed plane wave and a perturbation over 
CD only. The undisturbed plane wave, when integrated, gives a Dirac delta function 
in the forward direction. Thus the differential cross section comes from the sources 
on CD. The observation points are very far away so that to get the amplitude at an 
angle 6 we just sum all directions parallel to 6. 

Let p and be the polar coordinates in the plane (P). The amplitudes of the 
virtual sources in the disk CD arc then 

-i 1 + cos 9 



^disk(/3) 



A 



1 — c 
27r 



i2kB. 



1 



^i2ky/R^-p2 (n-1) 



1 



(C.l) 



The resulting scattering amplitude of Equation (jB.ip . now refered as /sphcrc(^), 
is obtained after moving the origin of the phases to the plane (P). This is simply done 
through a multiplication by c"'^*"'^ 

/sphere (0) = C-'2fc^ d0 ^disk (p) P e"='' ™^ ^ «dp 

(C.2) 

^i/cp cos (f> sin 9 



k 1 





cos 9 



2tt 



1 



J2ky/R^-p^ (n-1) 



dppdf/), 



'0 JO 

where the term p cos (f> sin 9 represents the path length difference between a source in 
the disk and a source at disk center for an observer placed at an infinite distance from 
the disk. 

But the integral over the angle is a Bcsscl function 

r2TT 



Jo{u) 



27r Jo 

Thus Equation ()C.2p becomes 



/sphere (6*) ifc 



1 



pjQ{kpsm9) 



1 



J2ky/R?-p^ (n-1) 



dp. 



(C.3) 



(C.4) 



This is as far as we can go with usual functions. Note that the scattering function in 
the forward direction is 



/. 



sphere I, 



= 0) = ifc 



pdp, 



which is exactly Equation (jB.4p . 
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We can get a closed form formula if the sphere is replaced by a disk of radius R 
and of height 2R along z. Equation (jC.4[) simplifies into 



/disk (6*) = ik 



= ik 



1 + cos ( 
2 

1 + cos ( 



pJo{kp sinO) 



,i2(n-l)fcfl 



1 - e 

R 



'\2kR. (n-1) 



dp 



pJo{kp sm6)dp 



(C.5) 







2{n~- l)kR 



.l + cos6' ■ 2{n-i)kR , . 
= 1 e 2 (— 2ijsm 



since the integral can be simplified by the fact that 



R 



Ji{kRsml 



uJo(u)du — RJi{R). 



Then the differential cross section reads 



do; 
d^ 



kR- 



1 



cos9 , . 2Ji(fci?sin( 

— sm [[n — l)kR\ 



(C.6) 



kR sin 6 

In order to compare predictions from Ramsauer approach and Mie, Equation (|C.6[1 
can be rewritten as a function of the size parameter x = kR and normalized by ttR^ 

2 



1 da 1 



7ri?2 d^ 



1 + cos 6* 



sin [x{n — 1)] 



2Ji(a;sin( 



X sm f 



(C.7) 



When the relative index of refraction n is equal to 1 , the differential cross section 
is null which is normal. Also, if the phase shift in the sphere is not far from a multiple 
of 27r, then again the cross section is null. This is the Ramsauer effect again: when 
the outside and the inside are in phase, they are producing an almost invisible sphere. 

Figure ICll shows the differential cross sections from Mie and from Ramsauer for 
a drop of water of 45.45 ^m radius, in air, at an incident wavelength of 0.6 /zm. As in 
the case of the total cross section, we find that the Ramsauer formula for a disk is a 
fair approximation of Mie for a sphere. 
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Figure CI. Angular scattering diagram for 45.45 fj,m droplet (n = 1.33) 
as obtained by Ramsauer (thick line) and Mie (dashed line) (calculated 
according |14) ) approaches for angles between 0° and 1.5°. The wavelength taken 
is equal to 0.6 ^m. 



